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ABSTRACT 



m 

We present a detailed re- analysis of the two ASCA Performance Verification observa- 
tions of the nearby Seyfert 1 galaxy MCG— 6— 30— 15. Confirming the results of Fabian 
et al. (1994), we find definite evidence for the O VII and O VIII K-shell absorption edges 
of the warm absorber and a doubling of the warm absorber column density within the 
3 weeks separating the two observations. No intra-day flux- correlated variability of the 
warm absorber is found. However, we report the discovery of an 'event' in which the 
warm absorber parameters temporarily change for ~ 10 000 s before returning to their 
original values. Possible interpretations are discussed but a contradiction remains: the 
constancy of the ionization state of the warm absorber argues that it lies at large dis- 
tances from the central source whereas the short term change in column density argues 
for small distances. Fluorescent iron emission is examined. As found by Fabian et al. 
(1994), the iron line is broad and strong (equivalent width ~300eV). The line profile 
is also suggestive of it being skewed. Such a line would be expected from a relativistic 
accretion disk. 

We also find very rapid primary X-ray variability. Assuming relativistic beaming to 
be unimportant, the derived efficiency is comparable to the maximum obtainable from 
accretion onto a Schwarzschild black hole. Correlated variability outside of the energy 
range of ASCA might exceed this maximum, thus requiring efficient accretion onto a 
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^ Key words: galaxies: individual: MCG— 6— 30— 15, galaxies: Seyfert, X-rays: galaxies, 

galaxies: active, plasmas, line: formation. 



1 INTRODUCTION 

The nearby bright Seyfert 1 galaxy MCG— 6— 30— 15 has 
been extensively studied in X-rays. A short EXOSAT ob- 
servation revealed X-ray variability on timescales of 1 500 s 
(Pounds, Turner & Warwick 1986), making it one of the 
first active galactic nuclei (AGN) in which such rapid vari- 
ability was known. Further EXOSAT observations showed 
evidence for the fluorescent K-shell emission and absorp- 
tion of cold iron (Nandra et al. 1989) expected from the 
illumination of dense, cold material located near the cen- 
tral engine by the primary X-rays (Guilbert & Rees 1988; 
Lightman & White 1988). Confirmation of these cold iron 
features by Ginga as well as the discovery of the associ- 
ated reflected continuum supported this picture (Nandra, 
Pounds & Stewart 1990; Matsuoka et al. 1990). Nandra, 
Pounds & Stewart (1990) also found evidence for partially- 
ionized, optically-thin gas along the line of sight from the 



iron K-edge. This material has become known as the warm 
absorber (Halpern 1984). Supporting evidence for the warm 
absorber model was provided by ROSAT Position Sensitive 
Proportional Counter (PSPC) observations which suggested 
the presence of absorption edges due to O VII and O vm at 
0.74 keV and 0.87 keV respectively (Nandra & Pounds 1992). 

ASCA (Tanaka et al. 1994) observed MCG-6-30-15 
twice during the Performance Verification (P V) phase of the 
mission. These observations are reported by Fabian et al. 
(1994), hereafter F94. They clearly resolved O VII and O vm 
K-shell edges for the first time. In addition, they found that 
the iron K-line at 6.4 keV was broad. This paper presents a 
detailed analysis of the ASCA PV data for MCG-6-30-15. 
Soft (< 2 keV) spectral features and variability are consid- 
ered within the framework of warm absorber models. Our 
examination of the hard spectral features concentrates on 
the fluorescent cold iron K-line. 
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Figure 1. SISO light curve (0.5-10 keV) for the July dataset with 
50 s bins. This light curve displays dramatic variations which are 
also seen in the other detectors. 



Section 2 briefly describes the observations and demon- 
strates the rapid variability exhibited by this object. Section 
3 examines the soft spectral features and variability within 
the framework of warm absorber models. Implications of the 
observations on the location and origin of the material are 
discussed. This is of wide importance: the study of these 
ionized absorbers is the first case in which a varying pho- 
toionizing continuum and the responding absorption edges 
are simultaneously and directly observed. Section 4 investi- 
gates the iron line feature, including its variability. Section 
5 provides a summary. 



2 OBSERVATIONS 

MCG-6-30-15 was observed by ASCA on 1993 July 8-9 
and 1993 July 31-August 1 (hereafter referred to as the July 
and August datasets). Each of these two observations extend 
over one day and have ~30 000 s of good data per detector 
(once contamination from the Earth and the South Atlantic 
Anomaly has been removed). Unless otherwise stated, Solid- 
state Imaging Spectrometer (SIS) data taken in both Faint 
and Bright mode have been combined. 

The SISO 0.5-10 keV light curve for the July dataset is 
shown in Fig. 1. The average 2-10 keV flux for this dataset 
was 4.6 x 10~ n erg cm~ 2 s _1 . The corresponding luminos- 
ity is 1.2 x 10 43 ergs _1 (assuming H = 50 km s" 1 Mpc" 1 ). 
Several rapid changes in luminosity can be seen. In the 
most dramatic event (Fig. 2), the 0.5-10 keV luminosity 



increases by 50 per cent (from ~ 1.8 x 10 ergs - to 
~ 2.8 x 10 43 ergs _1 ) in ~100s. The change in luminosity 
over the time taken for the change is 

^ 1n 41 -2 

— — « 10 ergs 



Figure 2. SISO light curve (0.5-10 keV) of the period between 
5 000 s and 10 000 s after the beginning of the observation. This 
shows one of the most dramatic changes in luminosity in which 
the SISO count rate increases by a factor of 1.5 in ~100 s. Similar 
behaviour is seen in the other detectors. As in Fig. 1, 50 s bins 
are used. 
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Figure 3. SISO light curve (0.5-10 keV) for the August dataset. 
50 s bins are used. Several strong flares are seen with the lumi- 
nosity doubling and then falling back to its original value within 
-10 000 s. 
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Assuming this is a true representation of the intrinsic 
variability (i.e. there are no effects from relativistic beam- 
ing), the efficiency of conversion of accreted mass into radia- 
tive energy (Fabian 1979) is at least 5 per cent. This is com- 
parable to the maximum efficiency obtainable from accretion 
onto a Schwarzschild black hole (Shapiro & Teukolsky 1983). 
This behaviour is also seen in the other SIS (SIS1) and the 
two Gas Imaging Spectrometers (GIS2 & GIS3). Correlated 
flux variations at energies outside of the ASCA band will 
serve to increase this lower limit to the efficiency. Thus, the 
efficiency may well exceed the maximum obtainable from a 
Schwarzschild black hole indicating the presence of highly 
efficient accretion onto a Kerr hole. 

Fig. 3 shows the SISO light curve for the August 
dataset. The average 2-10 keV flux for this dataset is 3.9 x 
10 -11 erg cm -2 s _1 . The corresponding 2-10keV luminosity 
is 1.1 x 10 43 ergs -1 . This lacks the dramatic short term vari- 
ation of the July dataset but does demonstrate strong 'flares' 
in which the luminosity doubles and then returns to the qui- 
escent value during a period of ~10 000 s. 



3 THE WARM ABSORBER 

F94 confirmed the presence of partially-ionized optically- 
thin material in the line of sight to the primary X-ray source. 
This has a large effect on the soft (0.5-2 keV) spectrum of 
MCG— 6— 30— 15, with absorption edges of O VII and Ovm 
being the dominant features (see Fig. 1 of F94). They also 
found evidence for a doubling of the column density of warm 
material between the two PV observations. Here we per- 
form a more detailed time resolved study of the warm ab- 
sorber. We examine variability within each dataset as well 
as changes between the datasets. 

3.1 Modelling the warm absorber 

Following F94, the photoionization code cloudy (Ferland 
1991) was used to examine the passage of the primary X- 
rays through an optically-thin, geometrically-thin shell of 
material, cloudy solves the equations of thermal and ion- 
ization equilibrium self-consistently with the radiation field. 
The incident spectrum is assumed to be a power law ex- 
tending between 13.6 eV and 40 keV with photon index V. 
We also constrain the density to be fixed at n = 10 9 ' 6 cm -3 , 
although the results are very insensitive to the actual density 
chosen (this is true for densities between n ~ 10 cm~ 3 and 
n ~ 10 11 cm -3 ). Such models can be parameterized by the 
column density, Nw, and ionization parameter*, £, of the 
material. These models are one-zone models in the sense 
that the density and ionization parameter of the absorbing 
gas is uniform. 

A three-parameter grid of such models was computed 

* The ionization parameter, £, is defined by 
? nB? 

where L is the incident ionizing luminosity, n is the electron num- 
ber density and R is the distance of the gas from the source of 
ionizing luminosity. This form for the ionization parameter is cho- 
sen due to the ease with which it can be related to observation. 



model 


July 


August 


parameter 


dataset 


dataset 


r 


1.98 ± 0.01 


1.87 ±0.01 


TOVII 


0.53 ± 0.04 


0.63 ±0.05 


T"OVIII 


0.19 ± 0.03 


0.44 ±0.04 


X 2 /dof 


679/470 


758/470 



Table 1. Results of power-law and Ovu/Ovm edge fit to the 
July and August datasets. Data from SISO and SIS1 were fit- 
ted simultaneously (with different instrumental normalisations) 
to obtain these results, r is the photon index of the primary 
power-law. tovii an d toviii are the optical depths at thresh- 
old of the O vn and O vm K-shell absorption edges respectively 
(with the edge energies held fixed at the physical values). The 
errors quoted are at the 90 per cent confidence level, Ax 2 = 2.7. 

in which Nw, £ and T, are free parameters. These models 
can be compared with data using the xspec atable utility. 
Such a procedure gives the best fit values for Nw, £ and T. 

A major limitation of these models is the restriction 
to conditions of thermal and ionization equilibrium. Indeed, 
the data presented here lead us to suggest that both ther- 
mal and ionization equilibria may not be satisfied. Krolik & 
Kriss (1995) have recently examined models of X-ray heated 
winds with the assumption of thermal equilibrium relaxed. 
They show that the relaxation of this assumption has an 
important influence on the resulting emission, reflection and 
absorption. However, the current equilibrium models do pro- 
vide a useful parameterisation of the data, provided they are 
interpreted carefully. 

3.2 Comparison of the July and August datasets 

In order to study spectral changes between the July and 
August datasets, a single integrated spectrum was extracted 
from each dataset for each of the four instruments on ASCA. 
The spectra were background subtracted using source free 
regions of the same observation. Initially, the data were fit- 
ted with a simple model consisting of a power law primary 
spectrum (photon index T), cold absorption (fixed at the 
Galactic value of N H « 4.1 x 10 20 cm- 2 ; Elvis, Wilkes & 
Lockman 1989) and the K-shell absorption edges of O vn 
and O vm (with the edge energies held fixed at the physical 
values of 0.74 keV and 0.87 keV respectively). [Throughout 
this work, energies of spectral features (e.g. absorption edges 
and emission lines) are quoted in the source rest-frame (with 
z = 0.008).] The free-parameters in such a fit are T, the opti- 
cal depths of the edges and the absolute normalisation. Data 
from the two SIS detectors were fitted simultaneously. The 
results from these fits are presented in Table 1. It can be 
seen that both edges undergo an increase in optical depth, 
although the change is more dramatic for the O vm edge. V 
also decreases by ~ 0.1. 

Self-consistent warm absorber models (as described in 
Section 3.1) were fitted simultaneously to all four instru- 
ments (leaving the normalisation of each instrument free 
to allow for the slight discrepancy between the four instru- 
ments). Fig. 4 shows the July dataset fitted with such a 
model. When performing such fits, cold absorption is in- 
cluded with a column density fixed at the Galactic value. 
Table 2 gives the results of such fits to the two datasets. As 
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model 
parameter 


July 
(full spectrum) 


August 
(full spectrum) 


July 

(< 4keV) 


August 
(< 4keV) 


L abs (0.5-2keV) 


0.8 


0.5 


0.8 


0.5 


L abs (2-10keV) 


1.2 


1.1 


1.1 


0.9 


£ un abs(0.5-2keV) 


1.1 


0.8 


1.2 


0.9 


i U nabs(2-10keV) 


1.2 


1.1 


1.1 


1.0 


log N w 

r 

X 2 /dof 


21.78±0.03 

32±2 
1.97±0.02 
1732/1390 


22.07±0.03 

43t^ 
1.87±0.02 
1637/1330 


21 7 n+o.02 

zl - ' -0.03 

26±2 
2.06±0.02 
926/825 


oo m+0.02 
zz.u3_ 03 

34±3 

2.00±0.02 

822/766 



Table 2. Results of three-parameter warm absorber fits to integrated July and August datasets. L a b s is the observed luminosity (after 
reprocessing by the warm absorber) whereas £ una bs is the inferred luminosity prior to warm absorption. All luminosities are given in 
units of 10 43 ergs -1 . T^w is given in units of cm -2 and £ is in units of erg cms -1 . The errors are quoted at the 90 per cent confidence 
level, Ax 2 = 2.7. Fits which ignore data above 4keV are free of any possible biasing introduced by the iron line. 
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Figure 4. Integrated (background subtracted) spectrum for the 
July dataset. All four instruments are used in the spectral fitting 
but, for clarity of presentation, only SIS0 (plain crosses) and GIS2 
(crosses with squares) are shown in the figure. Also shown is the 
best fit warm absorber model. 



can be seen from Fig. 4, the major discrepancies between 
these warm absorber models and the data occur in the range 
5-7 keV and may be interpreted in terms of K-shell emission 
of iron. Details of this feature are to be investigated in Sec- 
tion 4. It is relevant to the warm absorber analysis only in 
so far as it could bias the fit towards lower values of I\ To 
eliminate any such possible biasing, we performed a second 
fit including only these data below 4 keV (also shown in Ta- 
ble 2). Excluding data above 4keV leads to a steepening of 
the inferred photon index and a slight decrease in £ for the 
warm absorbing material. 

Comparing either the two full spectra or the two trun- 
cated spectra leads to the conclusion that the warm col- 
umn density has increased from (6.2 ± 0.4) x 10 21 cm~ 2 to 



Figure 5. Inferred £ plotted against the 2-10 keV luminosity for 
the July dataset. % and T are held constant at 6 x 10 21 cm -2 
and 2, respectively. Errors are given at the 90 per cent confidence 
level for one interesting parameter (Ax 2 = 2.7). 



(1.1 ± 0.1) x 10 22 cm~ 2 , i.e. the column density has doubled 
in three weeks. There is also evidence for an increase in £ 
between the July and August datasets despite a significant 
decrease in the ionizing luminosity. These results confirm 
the findings of F94. Comparison of the full spectra shows a 
decrease of T by 0.1. However, this effect is much less signif- 
icant in the truncated spectra and thus could be due to the 
behaviour of the hard spectrum. Some implications of these 
results are discussed presently. 



3.3 Short term variability 

The reaction of the warm absorbing material to the rapid 
fluctuations in the primary ionizing luminosity is an impor- 
tant diagnostic of the physical conditions within this gas. 
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model 
parameter 


July 
(first half) 


July 
(second half) 


August 
(first half) 


August 
(second half) 


L abs (0.5-2keV) 


0.9 


0.7 


0.5 


0.6 


iabs(2-10keV) 


1.3 


1.0 


0.9 


1.0 


£ U nabs(0.5-2keV) 


1.3 


1.0 


0.9 


1.0 


iunabs(2-10keV) 


1.4 


1.0 


1.0 


1.2 


log7V w 

r 

X 2 /dof 


21.75±0.03 

23 ± 2 

2 09+ 04 
z - u -0.03 

735/765 


21.89±0.06 

2.06±0.04 
583/570 


22.01 ± 0.03 

33 ±2 
2.02 ± 0.03 
674/675 


22 04 +0 03 
— 0.04 

35±2 
2.00 ±0.03 
681/680 



Table 3. Results of three-parameter warm absorber fits to each half of the July and August datasets. £ a b s is the observed luminosity 
(after reprocessing by the warm absorber) whereas L una bs is the inferred luminosity prior to warm absorption. All luminosities are given 
in units of 10 43 ergs -1 . -/V H is given in units of cm~ 2 and § is given in units of ergems -1 . Spectra for all four instruments were fitted 
simultaneously. Data above 4keV were ignored to prevent contamination by iron emission. The errors are quoted at the 90 per cent 
confidence level, Ax 2 = 2.7. 
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Figure 6. Summed SIS0 and SIS1 light curve in various bands 
for the July dataset. The panels show (from top to bottom) a) 
the 0.4-0.6 keV count rate, b) the 0.7-1.3 keV count rate, c) 2.5- 
5keV count rate, d) ratio of 2.5-5 keV data to 0.4-0.6 keV data 
(hardness ratio TZ\) and e) ratio of 2.5-5 keV data to 0.7-1.3 keV 
data (hardness ratio l^i)- All count rates are in counts a -1 . 250s 
bins were used. Note that during the period between 40 000 s and 
50 000 s H2 suffers a significant rise with no associated strong 
change in 1Z\. This is evidence for a rapid change in the warm 
absorber. 



Limits can be set on the photoionization and recombination 
timescales of the material thereby allowing the density to be 
constrained. Constraints on the density together with mea- 
sured values for £ and N-w give information on the distance 
of the material from the central ionizing source (Reynolds & 
Fabian 1994). In addition, inhomogeneities can be studied. 

Here we examine each dataset for short term variation 
of the warm absorber. Two techniques were used to analyse 
the short term variability of the warm absorber; variation of 



colour ratios and direct fitting of warm absorber models to 
spectra taken from various times within the dataset. 



3.3.1 The July dataset 

The July dataset was divided into two equal halves. Spec- 
tra were extracted for each half and warm absorber models 
fitted to the data below 4keV. Table 3 shows the results of 
such fits. There is evidence for an increase in both the col- 
umn density and ionization parameter of the warm material 
between the two halves of the dataset, despite a significant 
decrease in ionizing flux. The power-law index of the pri- 
mary spectrum, T, does not change significantly. 

The dataset was divided further. Spectra were extracted 
for each spacecraft orbit (giving 14 spectra for each instru- 
ment). Warm absorber models were then fitted to each of 
these spectra. Unfortunately, there were insufficient counts 
to usefully constrain three-parameter warm absorber mod- 
els. The major statistical degeneracy is a trade-off between 
Nw and £. To counter this degeneracy, we repeated the 
warm absorber fits holding Ny/ and T fixed at 6 x 10 21 cm -2 
and 2, respectively. The instrument normalisations were also 
tied together in an attempt to reduce statistical degenera- 
cies. Thus, £ and the overall model normalisation were the 
only free parameters in such fits. Fig. 5 shows the results 
of such fits. Here, we plot the inferred £ as a function of 
the 2-10 keV luminosity (taken to be indicative of the total 
ionizing luminosity). Within this parameterisation, signifi- 
cant changes of £ are seen. However, they are uncorrelated 
with ionizing luminosity and may be unrelated to the warm 
absorber. We can rule out the hypothesis of constant Nw, 
constant T and £ proportional to the ionizing flux. 

Changing colour ratios are another indicator of vary- 
ing warm absorbers (Netzer, Turner & George 1994). Warm 
absorption predominately affects the flux in the 0.7-1.3 keV 
range. On the other hand, 2.5-5 keV flux and 0.4-0.6 keV 
flux should be largely unaffected by warm absorption or iron 
emission. Thus, we define two hardness ratios 
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Figure 7. Confidence contours for the optical depths at threshold of the O vn and O vm edges (tovii and toviii respectively) during 
spectral fitting with the power-law and 2 edge model, (a), (b), (c) and (d) show orbits 4/5, 6/7, 8/9 and 10+ respectively. The contours 
show the 68 per cent, 90 per cent and 99 per cent confidence regions for two interesting parameters. 



For a fixed primary spectral shape, changes in the parame- 
ters of the warm absorber will lead to a change in 7?. 2 whilst 
leaving IZi unchanged. Fig. 6 shows the 0.4-0.6 keV count 
rate, the 0.7-13 keV count rate, the 2.5-5 keV count rate and 
these two hardness ratios. 

Fig 6 shows that for a period between 40 000 s and 
50 000 s after the beginning of the observation (i.e. space- 
craft orbits 8 & 9 of the observation), IZ2 undergoes a sig- 
nificant increase with little corresponding change in IZi . This 
'hardening event' was investigated further. SIS spectra were 
accumulated for the data taken during this period. Compar- 
ison spectra were accumulated for two periods before this 
event (orbits 4/5 and 6/7) and one period after the event 
(orbits 10+). Table 4 shows the results of various spectral 
fits to these spectra. Fig. 7 shows the confidence contours 
for the optical depths at threshold of the O VII and O vm 
edges (tovii and toviii respectively) in the power-law and 
2 edge fit. It is clear that the warm absorption undergoes a 
change (leading predominantly to an increase in toviii) and 
then returns to its original state. 
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3.3.2 The August dataset 

As with the July dataset, the August dataset was divided 
into two equal halves, spectra were extracted for each half 
and warm absorber models were fitted to each half of the 
dataset. Table 3 reports the results of these fits. There 
was no significant change in any of the warm absorber pa- 
rameters between the two halves of this dataset (although 
the average luminosity does not change dramatically be- 
tween these two halves of the dataset). As with the July 
dataset, spectra for individual spacecraft orbits do not con- 
tain enough counts to constrain three-parameter warm ab- 
sorber models. Fixing Nw and T at 1.2 x 10 22 cm~ 2 and 1.9 
leads to no significant variation in £ between the fits. 

Fig. 8 shows the 0.4-0.6 keV count rate, the 0.7-1.3 keV 
count rate, the 2.5-5 keV count rate and the hardness ra- 
tios 72.1 and 7J.2- Both hardness ratios are consistent with 
being constant over the period of the August observation 
(constant model gives x 2 /dof=75/74 and 81/74 for TZi and 
IZ2 respectively). Thus the data are consistent with the hy- 
pothesis that the warm absorber remains unchanged over 
the duration of the August dataset. 



Figure 8. Summed SIS0 and SIS1 light curve in various bands 
for the August dataset. The panels show (from top to bottom) a) 
the 0.4-0.6 keV count rate, b) the 0.7-1.3 keV count rate, c) 2.5- 
5keV count rate, d) ratio of 2.5-5 keV data to 0.4-0.6 keV data 
(hardness ratio TZi) and e) ratio of 2.5-5 keV data to 0.7-1.3 keV 
data (hardness ratio TZi). All count rates are in counts s -1 . 250s 
bins were used. No significant changes in hardness ratios are seen. 

3.4 Discussion 

3-4-1 Summary of the observational results 

As reported by F94, there is a clear change in the character- 
istics of the warm absorption between the July and August 
datasets. The dominant change is an increase in toviii, al- 
though tovii may also have increased. This accompanies a 
significant decrease in ionizing luminosity. One-zone warm 
absorber models interpret this as a doubling of Nw and a 
small increase in £. 

The July dataset shows evidence for short term warm 
absorber variability. In an event lasting ~ 10 000 s, a slight 
decrease in ionizing luminosity is accompanied by a sig- 
nificant increase in toviii- Again, one-zone warm absorber 
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model Orbits Orbits Orbits Orbits 

parameter 4 & 5 6 & 7 8 & 9 10+ 



£abs (0-5-2 keV) 0.9 1.0 0.6 0.7 

£ a bs(2-10keV) 1.4 1.5 1.3 1.0 

£ U nabs(0.5-2keV) 1.1 1.4 1.0 0.9 

£ U nab s (2-10keV) 1.4 1.5 1.4 1.1 



T 2.02 ± 0.02 2.03 ± 0.02 1.94 ±0.04 2.03 ± 0.03 

tovii 0.49 ± 0.07 0.51 ± 0.09 0.60 ±0.14 0.58 ± 0.09 

toviii 0.20 ± 0.07 0.17 ± 0.08 0.32 ±0.11 0.17 ±0.07 

X 2 /dof 232/244 282/232 204/201 270/220 



\ogN w 21.76 ± 0.08 21.64+g;gg 21.86Ig "g 21. 71+^ 

£ 30 ±7 19±\ 28+1 21 ±4 

T 2.02 ± 0.03 2.08l° °5 i- 98 !^ 2.08 ± 0.04 

X 2 /dof 255/248 282/232 200/201 278/220 



Table 4. Results of spectral fitting to data taken from orbits 4/5, 6/7, 8/9 and 10+ of the July dataset. Two separate models were fitted 
to each period of data; a simple power-law plus Ovu/O vm edge model (as in Table 1) and three-parameter warm absorber models (as 
described in Section 3.1). Only SIS data below 4keV was used. 



models interpret this as a doubling of Nw- The warm ab- 
sorber parameters then return to their pre-event values. 
This change is very similar to the change seen between the 
July and August datasets suggesting a possible bi-modal 
behaviour. We note that there are similar decreases in lumi- 
nosity which are not accompanied by a change in the warm 
absorber. No significant change in the warm absorption is 
seen in the August dataset. 



3.4-2 Important timescales 

Initially, we make two assumptions; that the warm absorber 
is one-zone (in the sense of Section 3.1) and that the ioniza- 
tion state of the material is dominated by photoionization by 
the primary radiation field. Relaxation of these assumptions 
will be discussed below. Given these assumptions, there are 
at least three characteristic timescales relevant to a discus- 
sion of warm absorber variability; the variability timescale of 
the primary ionizing flux (t V ar), the recombination timescale 
of the dominant ions within the warm material (t r ec) and the 
photoionization timescale of the dominant ions within the 
material (t p h). The recombination timescale (to all atomic 
levels) for O vm is given approximately by 

t rec « 200ng 1 T°- 6 s 

where n = 10 9 ngcm~ 3 is the number density of the warm 
material and T = 10 6 TsK is the temperature of the warm 
material (Allen 1973; Shull & van Steenberg 1982). Using 
the definition of £, this can be expressed as 

tree « 200£ 2 .Ri6.L43 1 Tb- 6 S 

where £ = 10 2 ^2 erg cms -1 , R = 10 16 iii6 cm is the distance 
of the material from an source with ionizing luminosity L = 
10 43 L 43 ergs- 1 . 

We can also provide a simple estimate of t p h ■ If we as- 
sume the ionizing continuum has a power-law form with pho- 
ton index T = 2 extending between B m i n and that 
the K-shell photoionization cross-section is proportional to 
v~ 3 (above the threshold energy), the K-shell photoioniza- 



tion timescale is given by 

t _ 16nAR 2 E th 

0~thL 

where <r t h is the photoionization cross-section at the thresh- 
old energy E t h and 

A = "'(fe) 

For O vil, this evaluates to 

tph » 20 R1 6 Lm s 

Thus, t ph < t^ provided £ 2 T 6 °- 6 > 0.1. 

If tvar is longer than the other relevant timescales, the 
material would achieve photoionization equilibrium. The ob- 
servationally inferred £ would be proportion to the ionizing 
flux. The present data rule out this scenario. We note that 
in at least one Ginga observation the inferred ionization pa- 
rameter does vary in proportional to the flux at 1 keV (see 
Fig. 7 of Nandra, Pounds & Stewart 1990). However, these 
Ginga fits may be sensitive to other forms of spectral vari- 
ability such as changes in T or changes in the amount of 
reflection. 

More generally, the ionization parameter £ would be 
expected to respond to increases in the ionizing flux on the 
photoionization timescale, whereas £ would respond to de- 
creases in the ionizing flux on the recombination timescale. 
The observed ionization state would then depend on a 
weighted average of the past flux history. 

3-4-3 Implications for warm absorber: beyond the 
one-zone model 

The present data show no apparent flux- correlated changes 
in £ over a period of a day but an increase in £ over the three 
week period. [Note that the apparent discrepancy between 
the change in £ and the change in the ionizing flux (i.e. £ 
shows an increase whereas the ionizing flux decreases) can 
be reconciled if the ionizing flux flared up during the un- 
observed period.] Thus, we can conclude that both t rec and 
tph are longer than a day, but that t p h is shorter than (or of 
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the order of) 3 weeks. These constraints on the timescales 
can be converted into constraints on the distance of the ma- 
terial from the ionizing source. Thus, on the basis of the 
above one-zone models, we would place the material between 
~ 6x 10 17 L43 cm and ~ 3x 10 18 L43 cm. This is thought to be 
somewhat outside the broad line region (BLR). However, we 
do observe an event in which the warm absorber parameters 
undergo a change in 10 000 s or less. Within the one-zone 
models, this would argue that t p h < 10 000 s and thus that 
the material lies within ~ 2 x 10 17 L43 cm. Therefore, the 
one-zone models produce a clear contradiction. 

Realistic models for this material may require the con- 
sideration of inhomogeneities or stratification within the ma- 
terial, the possibility that the ionization state is not dom- 
inated by photoionization and the possibility that the ma- 
terial is far from equilibrium. We now qualitatively discuss 
these some of these complications. 

Inhomogeneities in the warm material would have an 
important effect on variability. Such inhomogeneities could 
arise from multiphase behaviour of the material leading to 
the formation of warm clouds. Multiphase behaviour would 
be particularly relevant for models in which the warm mate- 
rial was considered as a component of the emission line re- 
gions (see the discussion in Reynolds and Fabian 1994). Sup- 
pose the 10 000 s 'hardening event' seen in the July dataset 
corresponds to the passage of a warm cloud in front of the 
primary X-ray source. Light travel arguments together with 
the observed 100 s primary variability of this object imply a 
size of ~ 10 13 cm for the X-ray emission region. To traverse 
this region in 10 000 s, the warm cloud must have a velocity 
greater than ~ 10 9 cms _1 . This places the warm material 
coincident with, or within, the BLR (i.e. ~ 10 16 cm) pro- 
vided gravitational forces dominate. 

Until now, we have been treating the primary X-ray 
source as a varying point source. Realistically, the X-ray 
source will have a finite extension. For example, primary X- 
rays may originate from the inner accretion disk (and hence 
have an extent of 10 13 -10 14 cm) with isolated 'flares' pro- 
ducing the variability on timescales down to ~ 100 s. If the 
warm absorber is inhomogeneous on length scales of 10 13 - 
10 14 cm, the properties of the absorbing gas along the line 
of sight to different flaring region could be different. This 
might lead to complicated and chaotic changes in the fit- 
ted warm absorber parameters. If the X-ray emission were 
dominated by only two flaring regions, the bi-modality of 
the warm absorber seen in the present A SCA datasets could 
have a natural explanation since two lines of sight are in- 
volved. Within this scenario, recombination and photoion- 
ization timescales would no longer dominate a discussion of 
inferred warm absorber variability. 

Alternatively, stratification could be important. As an 
illustration of a simple stratified model, consider the situa- 
tion in which the warm absorber consisted of two spatially 
distinct regions; a hot inner region responsible for some por- 
tion of the O VIII absorption edge and a cooler outer region 
responsible for the O VII absorption edge and some part of 
the O VIII absorption edge. Suppose the outer region has 
long photoionization and recombination timescales (longer 
than 3 weeks) whereas the inner region has a recombination 
timescale shorter than three weeks. The small drop in flux 
between the two PV observations could cause the hot, inner 
region to partially recombine, therefore causing the optical 



depth of the O vm edge to increase without an associated 
change in the O VII edge (which is compatible with the ob- 
served changes noted in Section 3.2). When fit with a one- 
zone warm absorber model, we would infer an increase in £ 
and Nw (since the Ovm/Ovil ratio would have increased 
and the total absorption would have increased). However, 
this model has difficulty in explaining the short term vari- 
ability of the warm absorber unless we invoke an inhomo- 
geneous hot region. This example shows the one-zone warm 
absorber parameterisation could lead to misleading results. 

Finally, we discuss the possibility that photoionization 
does not dominate the ionization state of the material. The 
presence of an additional heat source for the warm material 
would lead to an equilibrium state lying to the left of the 
photoionization equilibrium curve on Fig. 3 of Reynolds & 
Fabian (1995). Mechanical heating driven by the primary ra- 
diation pressure is one such possible heat source. In this case, 
the ionization state of the material would be determined 
mainly by collisional effects. The response of the ionization 
state and column density of the warm material to variabil- 
ity in the primary flux would depend on the response of the 
mechanical heating to the flux variability and not directly 
on the primary ionizing flux. Optical depth changes in the 
absorber might lead to feedback processes operating on the 
mechanical heating. Thus the variability of the state of the 
absorber would be complicated, including the possibility of 
hystersis effects and instabilities in the absorbing material. 
We note that the opposite case, in which heating exceeds 
radiative cooling, is examined by Krolik & Kriss (1995). 



4 THE IRON LINE 

MCG— 6— 30— 15 shows clear evidence for K-shell emission 
from iron, as discussed by F94. The reflection of X-rays from 
optically-thick, cold material (e.g. the accretion disk or the 
molecular torus) would result in the 6.4 keV fluorescent line 
of iron. In principle, the profile of this line is an important 
diagnostic of the conditions within the central engine. If this 
line originates from reflection by an accretion disk around a 
supermassive black hole, it is expected to have a character- 
istic asymmetric double peak structure resulting from the 
effects of Doppler shifting, relativistic beaming and gravi- 
tational reddening. Alternatively, if the line originates from 
reflection by the molecular torus that is proposed to exist 
in the standard model of Seyfert galaxies, we would expect 
a narrow line due to the small velocity dispersion of the 
torus material. The variability of this emission in response 
to primary flux variations also gives important information 
on the location of the line emitting material. The line flux 
variations would expected to follow the primary X-ray flux 
variations but would be lagged and smeared by the charac- 
teristic light travel time of the line emitting region. 

Here we examine the profile and variability of this emis- 
sion feature. All spectral fitting in this section includes back- 
ground subtraction (with background taken from source free 
regions within the same observation). 

4.1 Spectral Analysis 

Spectral fits to the iron emission feature were found to be ex- 
tremely sensitive to the parameterisation of the underlying 
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Figure 9. Ratio of the July data to a model consisting of a power- 
law (including the effects of reflection) and two absorption edges 
(to model the warm absorber). The open squares are SISO data 
(Faint mode data only) and the filled dots are GIS2 data (note 
that the spectral fitting was performed using all four instruments 
simultaneously but we only show two detectors for clarity of pre- 
sentation). The iron emission feature can be clearly seen to be 
broad. 



continuum. To eliminate as many systematic errors as possi- 
ble, Bright mode SIS data was discarded for the purposes of 
detailed spectral fitting (i.e. only Faint mode SIS data and 
GIS data was used). This reduces the amount of good SIS 
data to ~20 000s for the July dataset and ~10 000s for the 
August dataset (per SIS detector). Hence the July dataset 
will provide better statistical constraints than the August 
dataset [Note that F94 use only Bright mode data and thus 
obtain better constraints from the August dataset] 

Initially, each dataset was fitted with a power law and 
two absorption edges at low energies (to model the warm 
absorber). The resultant fits give x 2 /dof of 1704/1658 and 
1386/1423 for the July and August datasets respectively. 
As noted in Section 3, the dominant residuals in such a 
fit are between 5 keV and 7 keV and are attributed to iron 
emission. Fig. 9 shows the ratio of the July data to this 
model. This suggests the presence of a broad skew emission 
line. The addition of a Gaussian emission line proves to be 
highly significant (with Ax 2 of 73 and 59 for the July and 
August datasets respectively.) For the July dataset, the line 
energy is E = 6.2 ± 0.2 keV, the line width is a = 0.7 ± 
0.2 keV and the equivalent width is EW=400 ± 100 eV. The 
August dataset gives E = 5.9±0.3keV, a = 0.8tot k eV and 
EW=500t?°° eV. 

Ginga spectra of MCG— 6— 30— 15 show spectral flatten- 
ing at high energies, interpreted as the presence of a reflected 
continuum superposed on the primary power-law. Incorpo- 
rating such a reflected continuum significantly affects the 
results. The datasets have insufficient counts to constrain 



Figure 10. 2.5-5 keV count rate, 5.5-6.5 keV count rate and the 
hard to soft ratio of these two count rates (H3) for the July 
dataset. SISO and SIS1 data have been summed to produce these 
light curves. 500 s bins are used. 



the parameters of the reflection (such as solid angle sub- 
tended by the reflecting material at the X-ray source or 
inclination of the reflecting material). Thus, we fixed the 
reflection to be face-on with the reflecting material sub- 
tending 27t steradians at the primary X-ray source (i.e. as 
expected if the reflection were from an accretion disk illumi- 
nated by an isotropic X-ray source lying above the plane of 
the disk). Fitting a reflected power-law (and two low energy 
edges to model the warm absorber) to each dataset results 
in x 2 /dof of 1673/1658 and 1358/1424 for the July and Au- 
gust datasets respectively. This improves the fit over the 
simple (unreflected) power-law by Ax 2 of 31 and 34 respec- 
tively. The addition of a Gaussian emission line again leads 
to a further significant improvement with A^ 2 of 36 and 
28 for the two datasets respectively. The line parameters 



for the July data are E = 6.2 ± 0.3 keV, a 



0.6 



+0.3 



keV 



and EW=260 ± 100 eV. The August dataset gives E 



5.9l°jkeV, a ■ 



0.7±o j keV and EW=300 ± 100 eV. 



4.2 Short Term Variability 

Due to the small count rate within this emission feature, the 
only method to examine short term variability of the iron 
emission is via colour ratios. On the basis of ratio plots such 
as Fig. 9, we would expect somewhat more than 20 per cent 
of the 5.5-6.5 keV counts to originate from the iron emission. 
Thus we define the hardness ratio 



n 3 



5.5 — 6.5 keV count rate 
2.5 — 5keV count rate 



(3) 



In particular, we seek to test the hypothesis that the line 
is constant during large changes in primary flux. Assuming 
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Figure 11. 5.5-6.5 keV count rate as a function of 2.5-5 keV count 
rate for the July dataset. SIS0 and SIS1 data have been summed. 
Each point represents 500 s of data. 



the fluorescence is driven by the primary flux, this would 
give information on the size of the emitting region. 



4-2.1 The July dataset 

Fig. 10 shows the 2.5-5 keV count rate, 5.5-6.5 keV count 
rate and 72.3, as a function of time for the July dataset. 
The hypothesis of a constant is unacceptable with 
X 2 /dof=46/21. Fig. 11 shows the same data with the 5.5- 
6.5 keV count rate shown as a function of the 2.5-5 count 
rate. We can constrain any constant component in the 5.5- 
6.5 keV count rate to be less than 0.001 counts s _1 using a 
constant plus linear parameterisation of this relation. This 
sets an upper limit to the equivalent width of any constant 
iron emission to be ~ 10 eV (at the 90 per cent confidence 
level, A% 2 = 2.7). However, this result weakens significantly 
if the point with the lowest 5.5-6.5 keV count rate is ex- 
cluded. The upper limit on the constant component then 
becomes 0.04 counts s _1 corresponding to an upper limit 
for the equivalent width of any constant iron emission to be 
~ 250 eV. 



4-2.2 The August dataset 

Fig. 12 shows the 2.5-5 keV count rate, 5.5-6.5 keV count 
rate and IZ3 as a function of time for the August dataset. A 
striking feature of these light curves is the constancy of the 
colour ratio during the significant flare which occurs 20 000 s 
into the observation. The ratio is consistent (at the 90 per 
cent level) with a constant (x 2 /dof=14/19; best fit ratio of 
0.097±0.008). Fig. 13 shows the 5.5-6.5 keV count rate as a 
function of the 2.5-5 keV count rate. We can constrain any 
constant component in the 5.5-6.5 keV count rate to be less 



Figure 12. 2.5-5 keV count rate, 5.5-6.5 keV count rate and the 
hard to soft ratio of these two count rates (K3) for the August 
dataset. SIS0 and SIS1 data have been summed to produce these 
light curves. 500 s bins are used. 



than 0.02 counts s _1 , corresponding to an upper limit of 
~ 300 eV to the equivalent width for any constant line (at 
the 90 per cent confidence level). 

4.3 Discussion 

Detailed investigation of the iron emission would require 
more counts in this feature than possessed by the present 
datasets. However, we can say the following. 

i) Without the inclusion of the reflected continuum, 
spectral fitting with a single Gaussian gives a line which 
is very broad (a ~ 0.7 keV) and strong (EW~400eV). Such 
a line would have to originate from rapidly moving mate- 
rial (with velocity dispersions of ~ 30, 000 km s -1 or more) 
which is overabundant in iron (George & Fabian 1991). The 
spectral flattening and iron K-edge introduced by the re- 
flected continuum mimics part of the broad emission line. 
Fixing the reflected continuum to be from a face-on reflec- 
tor subtending 2-7T steradians at the primary X-ray source 
(i.e. maximizing the reflected continuum from any disc-like 
geometry), the best fit line is also broad but is weaker 
(EW~250eV). 

ii) Fig. 9 is suggestive of an asymmetric line, as might 
be expected from a relativistic accretion disk around a su- 
permassive black hole. Additionally, the Gaussian line fits 
consistently produce a centroid energy less then 6.4 keV. 
This could result from fitting a Gaussian line to an intrin- 
sically asymmetric emission feature. Spectral fitting using a 
line profile of a relativistic disk does not give a significantly 
better fit. Thus, longer ASCA observations are required in 
order to collect the counts necessary to define the line pro- 
file. Since this work was originally submitted, the results of 
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Figure 13. 5.5-6.5 keV count rate as a function of 2.5-5 keV 
count rate for the August dataset. SISO and SIS1 data have been 
summed. Each point represents 500 s of data. 



the pre-event value. Within the context of one-zone pho- 
toionization models, the resulting upper limit on the pho- 
toionization timescale (and hence the upper limit on the dis- 
tance of the warm material from the central source) is incon- 
sistent with the lack of flux-correlated variability. We tenta- 
tively identify this event with the passage of a warm cloud 
across the primary X-ray source. Velocity arguments then 
place this cloud within the BLR. To reconcile this with the 
lack of observed flux-correlated variability, we would have to 
postulate that the material is far from equilibrium and/or 
that its ionization state is not controlled by photoioniza- 
tion by the primary X-rays. Spatial structure within the pri- 
mary X-ray source together with an inhomogeneous warm 
absorber might also lead to complicated variability of the 
observed warm absorber parameters. It is interesting that 
in this case, where the photoionizing and corresponding ab- 
sorption edges are observed simultaneously, agreement with 
a simple photoionizing model appears lacking. 

Fluorescent cold iron emission is examined. The cur- 
rent data suggest that the emission line is broad (corre- 
sponding to velocities of 30 000kms _1 or more) and strong 
(EW~ 250 eV or more). There is also evidence for asymme- 
try within the line. Fluorescent emission from a relativistic 
accretion disk is the most natural explanation of these facts. 
There are insufficient counts within this emission feature to 
usefully constrain line variability. 



such an observation have become available (Tanaka et al. 
1995) and confirm the current interpretation. 

iii) There are insufficient data to produce useful con- 
straints on the variability of the iron line. Again, longer ob- 
servations are required. 



5 SUMMARY 

MCG— 6— 30— 15 is a nearby bright Seyfert 1 galaxy with 
strong spectral features due to X-ray reprocessing. In par- 
ticular, Ginga shows it to have a strong reflection compo- 
nent (Nandra & Pounds 1994) and ASCA shows it to have a 
strong warm absorber with prominent O VII and O vm edges 
in the soft X-ray spectrum. We have analysed the data from 
two ASCA PV observations of MCG-6-30-15 for spectral 
features and variability associated with such reprocessing. 

Confirming the results of F94, we find a doubling of 
the warm absorbing column density from 6 x 10 21 cm -2 to 
1.1 x 10 22 cm -2 between the three weeks separating the two 
ASCA observations. We also find a small increase in the 
ionization parameter £ despite a significant decrease in the 
ionizing flux. 

We also examine short term variability of the warm ab- 
sorber. We find no flux- correlated variability over the dura- 
tion of each observation. Assuming the material is in pho- 
toionization equilibrium with the primary radiation, this im- 
plies that the material must lie between 3 x 10 18 cm and 
1.5 x 10 19 cm from the source of ionizing radiation. This 
places it outside the BLR. However, we report the discovery 
of short term warm absorber variability in the July dataset. 
In an event lasting ~ 10 000 s, toviii (and, to a lesser extent, 
tovii) undergoes a significant increase and then returns to 
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